Microscopic derivation of the pion coupling to heavy— light mesons 
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I. INTRODUCTION 



Quark models endowed with effective quark current- 
current microscopic interactions, consistent with the re- 
quirements of chiral symmetry, have provided, through 
the years, an invaluably useful tool towards the construc- 
tion of an effective low-energy theory for the strong inter- 
actions [1-5]. Besides regularizing the ultraviolet diver- 
gences of the theory, these effective theories bring in the 
necessary interaction scale needed to make contact with 
the hadronic phenomenology. An important feature of 
this class of models is the essential convergence of results 
and conclusions across a variety of possible forms for the 
confining kernel. They fulfill the well-known low energy 
theorems of Gell-Mann, Oakes, and Renner [2], Gold- 
berger and Treiman [3], the Weinberg theorem [6], and 
so on. Using this formalism it is also possible to give, 
for this class of models, an analytic proof of the low- 
energy theorems in the light-quark sector [7]. It turns 
out that the chiral angle — the solution to the mass-gap 
equation — remains the only nontrivial characteristic of 
such a class of models and it defines the latter completely. 

In the present paper, we derive the generalized Gold- 
berger-Treiman (GT) relation for the heavy-light sys- 
tems. The heavy-light mesons have received a lot of 
attention in the past few years due to the discovery of 
new narrow states in the cs family [8]. Nowak et al. [9] 
and Bardeen and Hill [10] have postulated that the spec- 
trum of such states should reflect the pattern of dynam- 
ical chiral symmetry breaking. Namely in the heavy- 
quark limit, the mass spectrum is expected to be deter- 
mined by the light quark and, were chiral symmetry to 
be exact, the heavy-light states of opposite parity would 
have become degenerate. It then follows that physical 
splitting between parity doublers could be related to the 
scale of spontaneous breaking of chiral symmetry, that is 
the quark condensate or, phenomenologically, to the con- 



stituent quark mass. As more states in the open charm 
sector have been reported, alternative pictures have been 
examined, making this sector a good testing ground for 
phenomenological models [11]. Recently, for example, a 
canonical, quark model description has been investigated 
and it was argued that the new states can in fact 
be described as quark model states [12]. Notwithstand- 
ing the phenomenological successes of the naive quark 
model, it suffers from the serious shortcoming of being 
unable to account for the physics of chiral symmetry 
breaking. It therefore fails, among other issues basically 
related to pion physics, to constrain the treatment of 
strong decays by the underlying bound state dynamics. 
This failure results, for example, in predictions for cou- 
plings to the ground-state pseudoscalars which do not 
satisfy PC AC or the Goldberger-Treiman relations [13]. 
Roughly speaking, quark models consistent with chiral 
symmetry can be thought of being evolutions of the naive 
quark model supplemented by the constraint of the mass- 
gap equation, and therefore they should provide the cor- 
rect framework to study such strong decays. Here we 
shall show that such type of models leads indeed to the 
GT relation for the pion coupling between the opposite- 
parity states. This is a non-trivial result, which can only 
follow from the simultaneous use of a chiral invariant in- 
teraction together with a consistent treatment of the pion 
and heavy-light meson dynamics. 



The paper is organized as follows. In the following 
section we recapitulate the chiral quark models and dis- 
cuss the application to the ground-state pseudoscalar 
and heavy-quark sectors. In Section III, we present the 
proof of the Goldberger-Treiman relation for the heavy- 
light mesons. Summary and a brief discussion are given 
in Section IV. Throughout the paper we concentrate on 
the = 0^ heavy-light states, but the results can easily 
be generalized to higher spins. 



II. CHIRAL QUARK MODELS 

In this chapter, we give a short introduction to the chi- 
ral quark model. The model is described by the Hamil 
tonian 

H = j d^xip^{x){-ia-S7 +ml3)tp{x) 



with the quark current-current ( J^(x) = ■i/i(x)7^^i/'('^)! 
interaction parameterized by the instantaneous confin- 
ing kernel K'^{x — y). An important feature of th( 
models of the class (1) is the remarkable robustness o: 
their predictions with respect to variations of the quarl 
kernel. Although quantitative results may vary for dif- 
ferent kernels, the qualitative picture described rcmainf 
essentially the same. This is specially true for those re- 
lations enforced by the mechanism of chiral symmetrj 
breaking that should be independent of the spatial de- 
tails of the confining kernel provided it brings a natura 
confinement scale (hereafter called A'o). The only re- 
quirement is that it should be chirally symmetric. Wt 
illustrate this nice feature by deriving the Goldberger- 
Treiman relation connecting the pion coupling to heavy- 
light mesons to the mass splitting between chiral dou- 
blets, for the simplest Lorentz structure of the inter- 
quark interaction in the Hamiltonian (1) compatible witl 
the requirements of chiral symmetry and confinemeni 
KfA^ - y) = ^"W^^oV^dx - y|) [1, 3, 5, 14], for ar 
arbitrary confining kernel V{r). We also note that anal- 
ogous Hamiltonian with linearly rising potential consid- 
ered in one time and one spatial dimension reproduces 
the 't Hooft model for 2D QCD [15] in the Coulomb (ax- 
ial) gauge (see, for example, the original paper [16] or the 
review paper [17] and references therein). 

A standard way of proceeding with the investigation of 
the model (1) is to consider the self-interaction of quarks 
separately, thus introducing the notion of the dressed 
quarks: 



where the quark amplitudes 
1 



(2) 



i(p) = [Vl + + (a • p)i/l - Sp] wo(p), 



K-P) = + sp - (a • p)^l - sp\ Wo(-P) 

(3) 

are parameterized with the help of the chiral angle (pp 
[1-3], Sp = sin(^p, Cp = costpp. It is convenient to define 
the chiral angle varying in the range — 7r/2 < tpp < 7r/2, 
with the appropriate boimdary conditions f{0) = ^, 
(fiip oo) for the physical vacuum. The equation 




FIG. 1: A typical profile of the chiral angle — solution to the 
mass-gap equation. Momentum p is measured in the units of 
the potential strength Kq. 



which defines the profile of the chiral angle — the mass- 
gap equation — follows from the requirement that the 
quadratic part of the normally ordered Hamiltonian (1) 
should be diagonal in terms of the dressed quark creation 
and annihilation operators [3]. The mass-gap equation 
then takes the form: 



mc„ 



■psp 



J (2^ 



:V{p-]s.)[SkCp-XCkSp], (4) 



where a; = (p • k). In Fig. 1, we give a typical profile of 
the chiral angle. 

The dressed quark dispersive law can be evaluated then 

as 



msj, 



cPk 
(27r)3 



^(p - k) [spSk + xCpCk] . 

(5) 

Wc turn now to bound states of dressed quarks to 
the quark-antiquark mesons. Each mesonic state in this 
model is described with a two-component wave function, 
and the bound state equation acquires the form of a sys- 
tem of two coupled equations. The details can be found 
in the works [2, 3], where the Bethe-Salpeter equation 
formalism was developed or in Ref. [18], where a sec- 
ond Bogoliubov-like transformation is performed over the 
Hamiltonian (1) in order to define explicitly the com- 
pound mesons creation and annihilation operators. For 
the purpose of the present research, it is sufficient to con- 
sider two particular cases of this general mesonic bound- 
state equation. 



A. The case of the chiral pion 

In the chiral limit, the bound-state equation for the 
pion at rest is known to be given by the mass-gap equa- 



3 



tion (4) with m = [3, 18]. In the random phase ap- 
proximation (RPA) both components of the pionic wave 
function coincide and are simply siniy9p [2, 3, 18]. Beyond 
the chiral limit, the pion state is given by (p^r = 0): 



r 



&t(k)- 



+ d(-k)- 



' Ir I, 



d\-k)X{k) 



6(k)F(k) 



(6) 



|0>, 



with |0) being the broken BCS vacuum. Nf and Nc are, 
respectively, the rmmber of flavors and colors. 
The normalized wave functions are given by 



y(k) 



2^^ 



(7) 



2iU 



[sk - m^Ak] 



with Afe satisfying 
2EkAk = Sk+Cp 



(2^ 



-V{\s.-f\)[SkSq+XCkCg\Aq. (8) 



B. The case of a heavy— light meson 

The case of a heavy light meson, with one quark being 
infinitely heavy, leads to considerable simplifications of 
the bound-state equation in view of the fact that the neg- 
ative energy-spin component of the meson wave function, 
responsible for the time backward motion of the quark- 
antiquark pair, which is non-negligible for the light-light 
system, goes to zero with the bare mass increase of just 
one quark of the pair (heavy light systems), vanishing 
in the infinite mass limit of that quark. The resulting 
equation takes the form of a Schrodinger-like equation 
[19] 



^y(p-k) [c^Ck+ 



{a ■ p)((T • k)5p5fe *(k) = £;*(p). 



(9) 



with 



^j> = ~;j^ [VI + - Vi-cpJ • 



(10) 



Here the wave function is written as a product of the 
orbital part and the spin part, ^'(fc) = //\/2 • ^'(fc) for 
the = 0- and ^'{k) = (cr • k)/V2 • ^'(fc) for the = 
0"*" heavy-light meson, respectively. The conventional 



normalization of the wave function (consistent with the 
relativistic normalization) is 



where, in the heavy-quark limit, M approximates the 
mass of the heavy quark, and the mass of the heavy- 
light meson is given \yy M -\- E. 



III. GOLDBERGER-TREIMAN RELATION 
FOR HEAVY-LIGHT MESONS 

In this chapter, we focus on the details of the pion in- 
teraction with heavy light mesons. Consider the hadro- 
nic process which involves the chiral pion and the two 
heavy-light mesons for which, for the sake of trans- 
parency, we consider the D{JP = 0+) and the D'(J^ = 
0~) mesons with the wave functions ^{k) and ^'{k), re- 
spectively. 



A. The macroscopic derivation of the 
Goldberger— Treiman relation in QCD 

Let us consider the transition D — > D'-k. Then, using 
the relation 



(0lA«(0)l7r''(q)) = ^^g^<5»^ 



(12) 



we introduce the nn'-K coupling gnn'-K through the rela- 
tion 

{n'\Al\n) = (nV;i|n)nonpion - ^^^^i^^'^^"^' 

(13) 

where the subscript "nonpion" denotes the contribution 
free of the pion pole. Here D' , D are unit isospin doublets 
representing the flavor of the heavy-light n' and n state, 
respectively. 

On the other hand, the full matrix element of the axial 
current, corresponding to the l.h.s. of Eq. (13) can be 
written as: 



{n'\Al\n) = [{Pl + P,)GA{q^) 

- {P'-P^)Gs{q'')]D'^D, 



(14) 



with P' = Pn', P = Pn q = P' — P- To leading order in 
the heavy mass, conservation of the axial current in the 
chiral limit then demands 



2M{E' - E)Ga - q^Gs = 0. 



(15) 



Therefore, from Eq. (13) one can see that as q^ ^ 
it is possible to have Ga{0) ^ 0, if Gs is identified with 
the pion pole contribution. 



lim Gs{q') 

g2— ►O 



2- 4:Mf„g„ 



(16) 
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and Ga remains finite in the limit 0. The relation 

(15) then reads: 

1 



-{E' — E)Ga = fngnn'TT 



(17) 



A 

0+ 



c(-k) 



B 



0- 0+ 



c(-k) 



and is the sought Goldbcrgcr Trciman relation for the 
pion coupled to heavy-light mesons. 



V{p - k) I I V{p - k^ 

i|t(p)?/,(l<-)^g(p) \ \g (-p) g(-p) / /g(p) 



7r(X) 



7r(r) 



We now turn to the main subject of this paper and 
demonstrate how the Goldberger-Treiman relation (17) 
emerges microscopically in the chiral quark model. 

We are interested in matrix elements of the type 
{n'\0\n) where, in particular, O is the isovector axial 
current operator A°(0) = '0(O)7p^75V'(O), whereas \n) 
and \n') are the true cigenstatcs of the QCD Hamilto- 
nian Hqcd which, after integrating out gluon degrees 
of freedom, we approximate by the Hamiltonian of the 
chiral quark model (1). Formally, the Hamiltonian (1) 
is represented in the Fock space of single hadrons built 
on top of the RPA vacuum. It is convenient to split the 
Hamiltonian in two parts. 





c 




D 




B. The microscopic derivation of the 




c(-k) c(-p) 




c{-k) c(-p) 


Goldberger-Treiman relation in the chiral qucirk 




1 — * — ^ 




1 — " 1 


model 


0+ 


V{p - k) 


0- 0+ 


y(p-k) 



QCD 



{Hqcd -V) + V = H + V^H, 



(18) 



where V describes interactions with pions and H is as- 
sumed to produce the bare spectrum of hadronic states 
except for mass shifts. The bare part H is then identi- 
fied with those parts of the Hamiltonian (1), which, in 
the dressed quark basis (2), do not create (or annihilate) 
single light quark antiquark pairs. The latter are put 
into V and lead to non-vanishing matrix elements Vnn'TT 
between the opposite-parity meson states with one ad- 
ditional pion. Bearing in mind the low-energy limit, we 
thus assume that the spectrum can be well approximated 
by the bare spectrum (in the large- A^c limit) and the only 
relevant residual interactions are those with pions. In the 
interaction picture, the matrix elements of interest can 
be written as 

{n'\0\n) = {h'\U{+oo, 0)OU{0, -oo)|n), (19) 

where the states on the r.h.s. marked with tildes are 
eigenstates of the bare, pion coupling free, Hamiltonian, 
H and U is the evolution operator, 

i-tf 



U{tf,ti) =Texp 



■ (20) 



To leading order in the pion emission/ absorption (that 
is, in V) we obtain 



{n'\0\n) = {n'\0\h) + 

m 

+ E 



{n'\0\m){m\V\n) 
En — Em + ie 

{h'\V\m){m\0\fi) 
En' — Em + ie 



(21) 



9(k)\ \9(-k) g(-p)/ X(P) 



AX) 



FIG. 2: The four diagrams contributing to {D'n\V\D). 



The intermediate states always contain an extra pion. 
Now, because we are interested in the case O = A^(0), 
that is, in matrix elements of the type (0|^^(0)|7r''(k)), 
we can, for the case of the heavy-light n — > n'n transi- 
tions, saturate J2m ^ follows, {N = n, or n') 

(22) 

For the pion interacting with heavy-light mesons, the 
coupling constant Qnn'Tv can be introduced then as 

(5'7r"|l/|5) = 2Mignn'^D'^T''D{2nfd^^\P' + p^ - P). 

(23) 

and therefore, in order to derive the Goldberger Treiman 
relation, we are to evaluate microscopically the matrix 
element on the l.h.s. of Eq. (23). 

Obviously, the interaction potential V, when written 
in the quark basis, is given by the four diagrams shown 
in Fig. 2. In these diagrams c stands for the charm 
antiquark. The black dots correspond to the appro- 
priate spinor vertices — see Ref. [3]. As an example, 
we depict in diagram A one such vertex: u^{p)u{\s.) = 

i[l/l + Spy/1 + Sk + {<T ■ p){(T ■ k)^l - SpVl - SfeJ. As 

it was discussed before, the pion Salpeter amplitude has 
two components: the positive energy-spin component X 
and the negative energy-spin component Y. The four 
diagrams of Fig. 2 can be readily evaluated using the 
rules of Ref. [3]. For the chiral pion at rest, the diagrams 
A and B can be shown explicitly to cancel against each 
other so that one is left only with the diagrams C and 
D. This should not come as a surprise, once, in the chiral 
limit, pions at rest decouple from quarks. We use the 
notations Cx and Dy for the amplitudes corresponding 
to the matrix elements of V shown in Fig. 2 (diagrams C 
and D) with the interaction attached to the heavy quark. 



5 



Thus we have with the ampUtudes Cx and Dy calculated explicitly to 

be 

(£>V|y|D) = 2M [C^ + D^] (27r)3j(3)(P' + p,r - P) 

(24) 

I 



= 27m vWjTm^ J p^)5(sf^*'wv(<,-k)x(k)*(k)|yT^vm^-.yTTl;^/T^, 

(25) 

The high momentum (UV) behavior of integrals is determined by the short distance behavior of the kernel and 
wave functions. For purely confining kernels wave functions are soft and integrals are converging without need for 
renormalization. The short distance behavior of the kernel depends on integrating out gluons. In the approximation 
that keeps transverse gluons out from the Fock space the self-consistent solution the Dyson equation for the kernel 
leads to UV Coulomb potential with a running coupling softer then ~ 1/ log((?^) and leads to finite UV integrals. This 
is discussed, for example in [5, 20]. In the IR all integrals are finite for long range confining potentials for color singlet 
matrix elements as discussed, for example in [1-3]. Adding the two amplitudes together and using the definition of 
9dd't!i to leading order in m■,r^ we get: 

U9nn'. = ^^ j ^^*'(q)nq-k)^4^*(k)[yr^yrT7^-xyrT^yr^. (26) 

The non-pion contribution to the axial charge is computed from the expectation value of the time component of 
the axial current between the heavy-light meson states: 

Si?/ w (^'> 

Finally, integrating the Schrodinger-like (9) for ^(k) with ^''(k)cfc and, consequently, a similar equation for ^''(k) 
with \l/(k)cfc, we find: 



-[e'-e)Ga = ^^ j ^^*'(q)y(q-k)[c,(v^^T7;y^T^+xy^^^/^^) 

) 

(29) 



Ck{x^l + Sqyjl + Sfe + VI - Sk^/l - Sg)]*(k) 



or, after simple algebraic transformations: 



'-iE'^E)GA-^^ J ^^f'(q)nq-k)^^*(k)[yr^yTT^-,xyTT^VT^]. (30) 



Comparing Eqs. (26) and (30) we arrive immediately 
at the Goldberger-Treiman relation (17). 



IV. DISCUSSION 

In this paper, we proved the validity of the Goldberger- 
Treiman relation (17) in the chiral quark model. We as- 
sumed that the model Hamiltonian represents the fourth- 
compenet of a four- vector in a particular frame, the one 



in which the heavy mesons are at rest and can verify 
relations between matrix elements involving fourth com- 
ponent of vectors. Several important conclusions should 
be drawn from the presented consideration. First of all, 
it should be noted that had we dropped the Y piece 
of the pion wave function (as in naive quark models) 
we would have been not only violating the Goldberger- 
Treiman relation by 50% as should be expected, but also 
the two diagrams, A and B (with the interaction cou- 
pled solely to the light quarks), would have now survived 



6 



to give a result ineompatible with the Ga contribution. 
This is yet another manifestation of the Goldstone na- 
ture of the chiral pion which finds its natural implica- 
tion in the chiral model used in this work. Secondly, the 
Goldberger-Treiman relation can only be realized by the 
Salpeter solutions ^ and ^' of Ref. [19], which in turn 
contain the physics of progressive effective restoration of 
chiral symmetry for higher and higher excitations. More- 
over, with the Goldberger Treiman relation (17), we are 
in a position to go even further and to conclude that 
the pion coupling to excited mesons decreases as the ex- 
citation number of the heavy-light mesons increases 
according to the Goldberger-Treiman relation (17) with 
AM = M' - M and ^ 1. This gives an 

n — *oo n — *oo 

explicit pattern of the Goldstone boson decoupling from 
excited hadrons. Such a scenario was discussed in a re- 
cent paper [21], where the decrease of the pion coupling 
to excited hadrons was anticipated as a consequence of 
the chiral pion decoupling from dressed quarks. In the 



present work, this result comes out naturally in the mi- 
croscopic derivation of the Goldberger-Treiman relation 
and, at the same time, an explicit expression for the 
DD'tt coupling constant is found. 



V. ACKNOWLEDGMENT 

Work of A.N. was supported by the Federal Agency 
for Atomic Energy of Russian Federation, by the Fed- 
eral Programme of the Russian Ministry of Industry, 
Science, and Technology No. 40.052.1.1.1112, by the 
Russian Governmental Agreement N 02.434. 11.7091, 
and by grants RFFI-05-02-04012-NNIOa, DFG-436 RUS 
113/820/0-l(R), NSh-843.2006.2. A.P.S. would like to 
acknowledge support of the US DOE via the contract DE- 
FG0287ER40365 and to thank the support of the CFIF 
of the Instituto Superior Tecnico where part of this work 
was completed. 



[1] A. Amcr, A. Le Yaouanc, L. Oliver, O. Pene, and J.-C. 
Raynal, Phys. Rev. Lett. 50, 87 (1983); A. Le Yaouanc, 
L. Oliver, O. Pene, and J.-C. Raynal, Phys. Lett. B 134, 
249 (1984); Phys. Rev. D 29, 1233 (1984). 

[2] A. Le Yaouanc, L. Oliver, S. Ono, O. Pene and J.-C. 
Raynal, Phys. Rev. D 31, 137 (1985). 

[3] P. Bicudo and J. E. Ribeiro, Phys. Rev. D 42, 1611 
(1990); ibid., 1635 (1990); P. Bicudo, Phys. Rev. Lett. 
72, 1600 (1994); P. Bicudo, Phys. Rev. C 60, 035209 

(1999) . 

[4] S. L. Adler and A. C. Davis, Nucl. Phys. B 244, 469 

(1984); Y. L. Kalinovsky, L. Kaschluhn, and V. N. Per- 
vushin, Phys. Lett. B 231, 288 (1989); P. Bicudo, J. 
E. Ribeiro, and J. Rodrigues, Phys. Rev. C 52, 2144 
(1995); R. Horvat, D. Kekez, D. Palle, and D. Klabucar, 
Z. Phys. C 68, 303 (1995); N. Brambilla and A. Vairo, 
Phys. Lett. B 407, 167 (1997); Yu. A. Simonov and J. 
A. Tjon, Phys. Rev. D 62, 014501 (2000); F. J. Llanes- 
Estrada and S. R. Cotanch, Phys. Rev. Lett. 84, 1102 

(2000) ; N. Ligterink, E. S. Swanson, Phys. Rev. C 69, 
025204 (2004); A. Szczepaniak, E. S. Swanson, C. R.Ji, 
and S. R. Cotanch, Phys. Rev. Lett. 76, 2011 (1996); 
A. P. Szczepaniak and E. S. Swanson, Phys. Rev. D 55, 
1578 (1997); A. P. Szczepaniak and E. S. Swanson, Phys. 
Lett. B 577, 61 (2003). 

[5] A. P. Szczepaniak, E. S. Swanson, Phys. Rev. D 65, 
025012 (2002). 

[6] P. Bicudo, S. Cotanch, F. Llanes-Estrada, P. Maris, J. E. 
Ribeiro, and A. Szczepaniak, Phys. Rev. D 65, 076008 

(2002). 

[7] P. Bicudo, Phys. Rev. C 67, 035201 (2003). 

[8] BABAR Collaboration, B. Aubort et al, Phys. Rev. Lett. 
90, 242001 (2003); CLEO Collaboration, D. Besson et 
al., Phys. Rev. D 68, 032002 (2003); BABAR Collabora- 
tion, B. Aubert et al, arXiv:hep-ex/0607082. 

[9] M. A. Nowak, M. Rho, and I. Zahed, Phys. Rev. D 48, 
4370 (1993); M.A. Nowak, M. Rho and I. Zahed, Acta 
Phys. Polon. B 35, 2377 (2004). 
[10] W. A. Bardeen and C. T. Hill, Phys. Rev. D 49, 409 



(1994); W. A. Bardeen, E. J. Eichten, and C. T. Hill, 
Phys.Rev. D 68, 054024 (2003). 
[11] T. Barnes, F. E. Close, and H. J. Lipkin, Phys. Rev. 
D 68, 054006 (2003); A. P. Szczepaniak, Phys. Lett. B 
567, 23 (2003); E. van Beveren and G. Rupp, Phys. Rev. 
Lett. 91, 012003 (2003); K. Terasaki, Phys. Rev. D 68, 
011501 (2003); T. E. Browder, S. Pakvasa, and A. A. 
Petrov, Phys. Lett. B 578, 365 (2004); L. Maiani, F. 
Piccinini, A. D. Polosa and V. Riquer, Phys. Rev. Lett. 
93, 212002 (2004); Y. Q. Chen and X. Q. Li, Phys. Rev. 
Lett. 93, 232001 (2004); D. S. Hwangand D.W. Kim, 
Phys. Lett. B 601, 137 (2004); L. Maiani, F. Piccinini, 
A. D. Polosa and V. Riquer,Phys. Rev. D 71, 014028 
(2005); U. Dmitrasinovic, Phys. Rev. Lett. 94, 162002 
(2005). 

[12] F. E. Close, C. E. Thomas, O. Lakhina, and E. S. Swan- 
son, arXiv:hep-ph/0608139; O. Lakhina and E. S. Swan- 
son, arXiv:hep-ph/0608011. 

[13] F. E. Close and E. S. Swanson, Phys. Rev. D 72, 094004 

(2005). 

[14] A. P. Szczepaniak and E. S. Swanson, Phys. Rev. D 55, 

3987 (1997). 
[15] G. 't Hooft, Nucl. Phys. B 72 (1974) 461. 
[16] I. Bars and M. B. Green, Phys. Rev. D 17, 537 (1978). 
[17] Yu. S. Kalashnikova and A. V. Ncfediev, Yad. Fiz. 62, 

359 (1999) [Phys. Atom. Nucl. 62, 323 (1999)]; Usp. Fiz. 

Nauk 172, 378 (2002) [Phys. Usp. 45, 347 (2002)]; Yu. S. 

Kalashnikova, A. V. Nefediev, A. V. Volodin, Yad. Fiz. 

63, 1710 (2000) [Phys. Atom. Nucl. 63, 1623 (2000)]. 
[18] A. V. Nefediev and J. E. F. T. Ribeiro, Phys. Rev. D 70, 

094020 (2004). 

[19] Yu. S. Kalashnikova, A. V. Nefediev, and J. E. F. T. 

Ribeiro, Phys. Rev. D 72, 034020 (2005). 
[20] A. P. Szczepaniak and P. Krupinski, Phys. Rev. D 66, 

096006 (2002). 

[21] L. Ya. Glozman and A. V. Nefediev, Phys. Rev. D 73, 
074018 (2006). 



